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ABSTRACT 

Context. The second short run (SRa02) of the CoRoT space mission for Asteroseismology was partly devoted to stars belonging to 
the Mon OB2 association. An intense monitoring has been performed on Plaskett's Star (HD 47 129) and the unprecedented quality of 
the light curve allows us to shed new light on this very massive, non-eclipsing binary system. 

Aims. We particularly aimed at detecting periodic variability which might be associated with pulsations or interactions between both 
components. We also searched for variations related to the orbital cycle which could help to constrain the inclination and the mor- 
phology of the binary system. 

Methods. An iterative Fourier-based prewhitening and a multiperiodic fitting procedure have been applied to analyse the time se- 
ries and extract the frequencies of variations from the CoRoT light curve. We describe the noise properties to tentatively define an 
appropriate significance criterion and, in consequence, to only point out the peaks at a certain significance level. We also detect the 
variations related to the orbital motion and study them by using the NIGHTFALL program. 

Results. The periodogram computed from Plaskett's Star CoRoT light curve mainly exhibits a majority of peaks at low frequencies. 
Among these peaks, we highlight a list of about 43 values, including notably two different sets of harmonic frequencies whose funda- 
mental peaks are located at about 0.07 and 0.82 d~' . The former represents the orbital frequency of the binary system whilst the latter 
could probably be associated with non-radial pulsations. The study of the 0.07 d~' variations reveals the presence of a hot spot most 
probably situated on the primary star and facing the secondary. 

Conclusions. The investigation of this unique dataset constitutes a further step in the understanding of Plaskett's Star. These results 
provide a first basis for future seismic modelling and put forward the probable existence of non-radial pulsations in Plaskett's Star. 
Moreover, the fit of the orbital variations confirms the problem, already mentioned in previous works, of the distance of this system. 
The existence of a hot region between both components renders the determination of the inclination ambiguous. 
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1. Introduction 

Plaskett's Star, or HD47129, has long been considered as one 
of the most massive binary systems in our Galaxy. For nine 
decades, this star has not stopped fascinating by its complex- 
ity. It has been the subject of several investigations in different 
wavelength domains but, despite these efforts, the system is still 

not fully understood. 

Plaskett's Star is a non-eclipsing binary ( Morr ison! 1 1978I) 
composed of two very massive and luminous O-type stars. The 
secondary component features broad and shallow stellar lines 
suggesting that this star rotates rapidly. Its projected rotational 
velocity has been estimated close to 300 kms -1 while that of 



the pr imary has been measured at about 75 km s~' dLinder et alJ 
2008). As a consequence of this high rotation speed, the sec- 
ondary star probably presents a temperature gradient between 
the poles and equator which could bias the determination of its 
spectral type but, most importantly, the secondary has a rota- 
tionally flattened wind. This configuration probably accounts for 
several proper ties of the wind intera ction zone. The study of the 
Ha region bv lWiggs & Giesl d 19921) and llinder et all (120081) as 
well as the analysis of iLinder et al. (1200a) in the X-rav domain 
have confirmed this assumption. 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany and Spain. 
** Senior Research Associate F.R.S.-FNRS 
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iBagnuolo et al.1 (1 19921) applied a tomographic technique to 
the International Ultraviolet Explorer (IUE) data to separate the 
contribution of the secondary star from the primary spectrum. 
Their investigation of the individual spectral components has 
provided the spectral types of 07.5 I and 06 I for the primary 
and secondary stars, respectively, and a mass ratio of MzlM\ - 
1.18 ± 0.12. Assuming an inclination of 71 ± 9° as estimated 
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bv lRudv & Herman! (1 19781 ) from polarimetry, they found masses 
eq ual to 42.5 Mp for the primary and 51.6 M Q for the secondary. 

iLinder et alj (120081) found, from high-resolution optical 
spectra, minimum masses of 45.4 and 47.3 M Q for the pri- 
mary and secondary, respectively, implying a mass ratio of about 
1.05 + 0.05. They also disenta ngled the spectra by using an al- 
gorithm based on the method of lGonzalez & Levatol (12006). The 
individual spectra indicated an 08 III/I+07.5 V/III binary sys- 
tem. In add ition, these authors use d the model atmosphere code 
CMFGEN dHihier& Milled! 19981) to derive the wind and the 
photospheric properties of both components. The major point 
of their study is the confirmation of a N overabundance and a 
C depletion of the primary star, providing additional proof for a 
binary system in a post Roche lobe overflow evolutionary stage 
where matter has been transferred from the primary to the sec- 
ondary star. 

The most disturbing point concerning Plaskett's Star is how- 
ever the discrepancy between the luminosi t y of bo th components 
and their dynamical masses. ILinder et al.l ^2008) explained that 
these stars have spectral types which are too late for their masses. 
A solution to this problem is to assume a larger distance of the 
star although this would imply that Plaskett's Star does not be- 
long to the Mon OB2 association. Despite the numerous inves- 
tigations undertaken to probe the physics of this binary system 
and its compone nts, Plaskett' s Star still hides a part of mystery. 

Even though iReese et all d2009l) have emphasized the diffi- 
culty to detect pulsation modes in rapidly rotating s tars, there are 
some examples of rapidly rotating O stars ( C Qph. iKambe et al.l 
[l997t HD 93521 iHowarth & Reidl fl993l lRauwetal.1 12008) 
where spectroscopic and photometric variability, probably re- 
lated to non-radial pulsation modes with periods of a few 
hours, have been identified. These rapid rotators have prop- 
erties quite reminiscent of those of the secondary component 
in Plaskett's Star, which also displays line profile variability 
dLinder et al.l l2008) although the existing spectroscopic data of 
the system are too scarce to characterize properly these varia- 
tions. Asteroseismology could therefore provide new insight into 
the properties of the components of Plaskett's Star. 

For this purpose but also to further constrain the binary sys- 
tem itself, Plaskett's Star has been chosen as one of the C oRoT 
(Convection, Rotation and planetary Transits, Bag lin et al . 2006; 
lAuvergne et al.ll2Q09h satellite targets in the Asteroseismology 
field. The unprecedented quality of the CoRoT light curve 
clearly allows us to search for variations linked to the orbital 
period of the system and to determine the possible existence of 
very low-amplitude variations due to, e.g., the presence of radial 
and non-radial pulsation modes. 

The present paper describes a complete and detailed analy- 
sis of the light curve of Plaskett's Star observed by the CoRoT 
satellite. We organize it as follows. In Section [2] we present the 
CoRoT data. Section [3] is devoted to the frequency analysis and 
to a study of the noise properties of the data. In Section [4] we 
discuss the inclination of the binary system and some other re- 
sults obtained from the analysis of the light curve. Finally, we 
set forth our conclusions in Section[5] 



2. The CoRoT data 

Plaskett's Star was observed during the second CoRoT short 
run SRa02 pointing towards the anticenter of the Galaxy. 
The observations were made between October 08, 2008 
{HJD = 2454748.485467) and November 12, 2008 (HJD = 
2 454 782.8 19956). Hereafter, we will express the date as HJD- 
2450000. The time series obtained is spread over AT = 



34.334489 days, with a sampling of one point every 32 s. This 
involves a frequency resolution of 1 / AT = 0.029 1 25 d~ 1 . 
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Fig. 1. Spectral window of the final version of Plaskett's Star 
light curve, observed during the second short run of the mission. 



The raw light curve of Plaskett's Star contains 92696 points. 
At first, we discarded flagged points potentially corrupted by the 
instrumental conditions (e.g., the changes of CCD masks or the 
generation of outliers). These observations represent about 4.2% 
of the datapoints of the CoRoT light curve of Plaskett's Star. In 
addition, we searched the CoRoT light curve for possible jumps 
or the presence of discontinuities due to a change of the CCD 
temperature and we corrected all of them. Secondly, we also 
discarded the flagged points associated with the environmental 
perturbation of the satellite mainly due to the regular passage 
through the South Atlantic Anomaly (S AA) and other Earth or- 
bit perturbations (lAuvergne et al.ll2009l) . Since this passage oc- 
curs twice in a sidereal day, the spectral window (see the inset in 
Fig- [U) presents a structure composed of a first peak (4.6% of the 
amplitude) close to 2.007 d -1 and a second one (1.4% of the am- 
plitude) close to 4.014 d -1 . Moreover, gaps due to the orbital pe- 
riod of the satellite (6184 s) produce other structures with peaks 
around 13.972 d~' and their harmonics (Fig.Q}. The percentage 
of lost datapoints because of environmental conditions is only 
9.6% for HD 47129. In addition, the spectral window exhibits a 
89.0% peak at / = 2699.764 d -1 . This corresponds to the sam- 
pling regularity having the highest frequency. Consequently, a 
pseudo-Nyquist frequency can be located at fy y = 1349.882 d -1 
(leading to a step of 32.003 s). We use the word "pseudo" to 
point out the fact that the aliasing is not pure. 

An interpolation of all flagged points of the light curve, to fill 
the gaps, cleans the spectral windows of all these peaks but this 
process introduces systematic effects which could generate erro- 
neous values of the frequencies, affecting the scientific results. 
Accordingly, we decided not to interpolate among the remaining 
points of the light curve. 

All stars observed in the CoRoT field of view present a 
global slope in their light cu rve, with probably an instrumental 
origin (lAuvergne et al.l l2009 attributed this drift to the ageing of 
the CCDs). Even though the light curve of Plaskett's Star does 
not seem to be affected by it, maybe because of its high flux vari- 
ability, we decided to remove this long-term trend. In order not 
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to bias the possible orbital or long-term variations present in the 
light curve, the strictly linear trend appears to be the best choice 
among the different possible trends. However, as we will show 
in Section [3] the low frequencies, close to 0, still appear in the 
semi-amplitude spectrum after removing the trend. 

Finally, the CoRoT light curve has also been converted to 
magnitude from the expression m - -2.5 log(f ) + C where F 
gives the CoRoT flux expressed in e~ s _1 and C represents a cal- 
ibration constant. We have estimated this constant at a value of 
23.09 ± 0.01 mag after having compared the CoR oT magnitude 
to that in the V-band quoted bv lLinder et al.l(l2008l) . We note that 
this conversion in magnitude does not change the decomposition 
in frequencies reported below. The final version of Plaskett's Star 
light curve (Fig. |2j is composed of 79896 points. 

3. Analysis of Plaskett's Star CoRoT light curve 

3.1. Looking for periodic structures 

The CoRoT light curve confirms that Plaskett's Star is not 
an eclipsing binary although we detect variations that could 
be linked to the orb ital motion of the binary system (P 01 -b = 
14.39625 days, see iLinder et al] 120081) . The light curve also 
presents large amplitude peak-to-peak variations on shorter time 
scales. In order to perform a complete investigation of the 
frequencies present in the CoRoT light curve, we applied a 
Fourier analysis based on the Heck, Ma nfroid & Mersch (here- 
after H MM) technique dHeck et al.ll985l revised by|Gpsset et al. 
1200 lb . We emph asize that this meth od is comparable to the 
Ferraz-Mello one dFerraz-Mellol 1 98 H) and is specially designed 
to handle time series with unevenly spaced data. Moreover, its 
mathematical expression for the power spectrum has the ad- 
vantage of cor recting some deficiencies of other methods such 
as the one of IScargla (Q~982). The Fourier technique of HMM, 
used here, is equivalent, at ea ch individual frequ ency, to a least- 
squares fit of a sine function dGosset et alJl200lh . In the present 
paper, we will actually use the semi-amplitude spectrum instead 
of the classical power spectrum, implying directly that the ordi- 
nates represent the amplitude term in front of the sine function. 

The semi-amplitude spectrum (Fig. [3]> clearly shows that 
most of the power is concentrated in the / < 6.0 d -1 frequency 
domain and that aliases, generated by the satellite orbital cy- 
cle, are present between 11.0 and 17.0 d _1 . Although we re- 
moved the linear trend from the dataset, there remains signal 
close to frequency 0.0 d _1 . The spectrum is dominated by a set 
of frequencies distributed according to a regular pattern includ- 
ing the fundamental frequency at about 0.82 d -1 and its possible 
harmonics at about 1.64 d -1 , 2.46 d" 1 , 3.28 d _1 and 4.10 d" 1 . 
The 2.007 d -1 aliases of the two first (main) frequencies are 
also visible at / = 2.83 d -1 and 3.65 d -1 . Besides these fre- 
quencies, the semi-amplitude spectrum confirms the presence of 
/ = 0.07 d"\ / = 0.14 d 1 and probably / = 0.21 d"\ i.e., the 
orbital frequency and its harmonics which constitute a second 
set of frequencies. Further outstanding peaks are also present at 
/ ~ 0.35 d _1 , / ~ 0.65 d- 1 and / ~ 0.95 d _1 . The first one 
having a wider peak that expected from the time basis, it could 
actually be a blend of several frequencies (at least two, or one 
plus an alias of / = 1.64 d~'). 

In order to analyse the time series in a more systematic 
way, the in-depth determination of the different frequencies was 
done in two steps. First, we proceeded to a classical iterative 
prewhitening of the signal, frequency by frequency. At each step, 
the semi-amplitude spectrum is computed and the frequency 
value is selected by the position of the highest peak in the spec- 



trum. The amplitude and the phase, corresponding to this fre- 
quency, are directly computed from the Fourier function. We 
however initiated the process by removing the set of 5 frequen- 
cies composed, notably, of the two highest peaks found in the 
periodogram (/ = 0.82 and / = 1.64 d" 1 ) and the frequen- 
cies belonging to this set. Since the frequency / = 0.82 d~' 
presents possible harmonics, we decided to perform the fit on all 
of them together, taking into account the existing correlation be- 
tween the power at the different frequencies. We thus designed 
a generalized periodogram giving at each frequency the power 
included in a frequency and its harmonics. The method is based 
on a part icular case of the ex tension of the HMM method pro- 
posed bv lGosset et al.l (|2001). The generalized periodogram ex- 
hibits the / = 0.82 d _1 peak characterized by the natural width, 
thus further confirming the fact that this sequence of peaks are 
indeed true harmonics. These frequencies were removed alto- 
gether from the data. Next, the semi-amplitude spectrum has 
been computed again in order to detect the new highest peak 
to remove. The analysed data are prewhitened for this frequency 
and we repeated this process until reaching the noise level of the 
data. Fig.|4]compares the semi-amplitude spectra of the data and 
of the data prewhitened for the first 20 frequencies, the first 40 
and the first 60 ones. 

Although efficient, the iterative removing of several frequen- 
cies, one by one, by using the basic HMM method, can not be 
the final procedure. Since the light curve presents gaps and the 
sampling is irregular, the height of a peak is dependent on the 
height of other peaks in the periodogram. Therefore, we need 
to fit all the listed frequencies together. For this purpose, we 
used the extension of the HMM Fourier method to the multifre- 
quency adju s tment , a high order Fourier method introduced by 
iGosset et al.l d200ll their Eq. A13 to A19). Because of an exces- 
sive computation time, we applied the multifrequency analysis to 
binned data to be able to deal with some amount of frequencies 
at the same time. Indeed, since no peak is clearly dominant above 
10 d _1 , we have binned the data with a step of one twentieth of 
a day, i.e., with a pseudo-Nyquist frequency of fy y = 10 d~\ 
reducing the CoRoT light curve into a set of 684 points. The 
multifrequency algorithm thus takes into account the mutual in- 
fluence of the different frequencies and investigates, inside the 
natural widths of the individual selected peaks, several positions, 
i.e., refined values for each of the frequencies in order to find the 
current best fit. 

The results obtained by both methods confirm the detection, 
in the semi-amplitude spectrum, of a first structure composed of 
a main frequency and six harmonics. The fundamental peak is 
situated around 0.823 d _1 and its harmonics around 1.646 d -1 , 
2.469 d"\ 3.292 d 1 , 4.1 15 d" 1 , 4.938 d" 1 and 5.761 d" 1 , respec- 
tively. Furthermore, we also clearly detect the peak correspond- 
ing to the orbital period of Plaskett's Star (/ ~ 0.069 d -1 ) as well 
as two harmonics of this frequency. We also emphasize the de- 
tection of two frequencies, actually unresolved, at / = 0.368 d -1 
and / = 0.399 d 1 . 

The set of frequencies, computed in that way, represents the 
final set (listed in Appendix). Actually, we stop the iterative mul- 
tifrequency procedure when the noise level is reached (statisti- 
cally). This critical level is particularly difficult to estimate and 
constitutes the main topic of the next section. 

3.2. The noise properties and significant frequencies 

From the analysis of the time series, illustrated in Figs. [3] and [4] 
we conclude that there is a clear excess of power at low fre- 
quencies. It represents either the reality of the underlying de- 
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Fig. 2. Full detrended CoRoT light curve of Plaskett's Star, containing 79896 points, observed over about 34.33 days and converted 
to magnitude. All gaps present in the data are due to the flagged points described in Sect. [2] 
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Fig. 3. Semi-amplitude spectrum of the light curve of Plaskett's Star computed by the HMM method from the unfiagged CoRoT data. 
We clearly see the aliases due to the orbital period of the satellite near 1 1-17 d -1 . The inset shows a zoom-in on the semi-amplitude 
spectrum in the low-frequency domain. 



terministic signal or the fact that the signal is partly made of red 
noise, i.e., is partly stochastic. Visible on the log-log plot (Fig. [5} 
of the periodogram of Plaskett's Star, this excess of power at 
low frequenc ies (log f < 0.6) can be , for example, described, as 
suggested by Stanishev et al. (2002)), by using a function of the 
form: 



P(f) = 



C 



l + {27irf)y 



(i) 



with y being related to the slope of the linear part and t be- 
ing an estimation of the mean duration of the dominant tran- 
sient structures in the light curve. A least square fit of the semi- 



amplitude spectrum in the low-frequency domain yields param- 
eters y - 2.3 and t = 0.12 d~'. Assuming an origin at least 
partly intrinsic to the star, several processes could explain the 
behaviour of the stochastic component of the signal. It could be 
either due to an onset of clumping at the ste llar photosphe re or, 
if we make the paral lelism with the works of Harvev ( 1985]) and 
Aigr ain et al.l (120041) in heliose ismology, to some kind of granu- 
lation. Indeed, Cantiello et al. (2009) suggested that the convec- 
tion zone induced by the iron opacity bump can have an impact 
on the stellar surface behavi our and thus cou l d be responsible for 
the existence of red noise. Belk acem et al. further sug- 

gested that this iron convection zone could generate stochasti- 
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Frequency (d ) 

Fig. 4. For comparison, the top panel shows the semi-amplitude 
spectrum before prewhitening, the subsequent panels represent 
the semi-amplitude spectrum after prewhitening 20, 40 and 60 
frequencies, respectively. Red lines exhibit the critical level at 
the significance level of 0.01 under a null hypothesis of white 
noise. 



cally excited modes in massive stars. Finally, other alternatives 
are also mentioned to explain this ori gin as the non-linearity (for 
more details, see e.g., Perdang 2009). 
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Fig. 5. Log-log plot of the semi-amplitude spectrum of Plaskett's 
Star. The red curve represents the fitted function (Eq.[TJ. 



Despite the possible presence of red noise, it is thus impor- 
tant to check whether peaks detected in the semi-amplitude spec- 
trum might result from random variations rather than represent- 
ing a periodic signal. However, the evaluation of the absolute 
significance of a peak in a periodogram of a time series with 
an uneven sampling is a c ontroversial problem (e.g jRauw et al.l 
2008: lFrescura et al.l2008h . Furthermore, the detection threshold 
level will be affected by this assumption of coloured noise, mak- 
ing its estimate more difficult. As a first order indication, we have 
thus applied a statistical test to establish whether the frequen- 
cies are significant, especially for those with smaller amplitudes. 
The probability that the power at several inspected frequencies 
exceeds a threshold z under the null hypothesis of a stochastic 



process of variance < r 2 f (func tion of frequencies) is given by the 
empirical formula of lGosse l J2007h : 



Prob[Z max > z] = 1 - (e 



where Z„ 



. e -0.93z+ln(0.8W ) 



(2) 



Pf/Vf with Pf the power in 



= max Z(/), j 

0</<Ay 

the power spectrum at the frequency / which is related to the 
semi-amplitude by Pf = NqA 2 JA and No represents the num- 
ber of datapoints. Although this expression is only adapted to 
even sampling, it remains a good approximation (since no bet- 
ter alternative exists) for uneven ones, especially for the CoRoT 
sampling which is more reminiscent of a regular but gapped one. 
We adopted the number of datapoints in the binned light curve, 
i.e., A^o = 684. A first estimation of the significance level can 
be made by adopting a null hypothesis of white noise (i.e., crj. 
not depending on /). In this case, we detect 79 frequencies (60 
frequencies) at a significance level of 0.01 (0.001) or lower. A 
second approach is to adopt a more complex stochastic process 
as null hypothesis. In this case, we defined a stochastic distribu- 
tion law by assuming an empirical function based on white noise 
above / = 6.3 d~' and on red noise defined by the continuum in 
the power spectrum below / = 6.3 d _I . We have normalized the 
above fitted function in order to define a properly scaled <xj func- 
tion. Under this refined null hypothesis of red noise, the number 
of significant frequencies decreases to about 43 (38) at a signif- 
icance level of 0.01 (0.001). Table Q] exhibits the details of the 
significant frequencies detected under the null hypothesis of red 
noise. The columns represent the sequence number of the fre- 
quency, the frequency itself, its amplitude, the significance level 
under the null-hypothesis of white and red noise and some com- 
ments, respectively. 

If we look at the evolution of the data variance with the itera- 
tive prewhitening process (Fig.[6]l, we note that there is no jump 
in the curve, which underlines the difficulty to define a thresh- 
old. Indeed, both structures formed by the fundamental peaks at 
/ = 0.823 and / = 0.069 d _1 already explain 72% of the vari- 
ance of the CoRoT light curve and, after 20 frequencies, the total 
variance of the prewhitened observations strongly decreased. 

The decomposition of the signal into frequencies is a formal 
process that does not necessarily represent the physical truth. 
Moreover, it neglects non linearities and complex phenomena. 
It is thus impossible to derive a list of frequencies free of con- 
taminations. However, the truly existing frequencies (orbit, pos- 
sible pulsations,...) should be part of this derived list (reported 
in the Appendix). We also decided to use an additional method 
to contribute to the determination of the reliable frequencies al- 
though this method is not perfect either. We split the data into 
two halves: a first one gathering the observations taken between 
HJD = 4748.48 and 4765.63 and a second one taking into ac- 
count the data collected between HJD = 4765.63 and 4782.82. 
We applied the frequency research on each sample separately by 
using the same technique as above and we directly compared 
both frequency lists. This method allowed us to put forward a 
list of 30 frequencies common to both dataset. 

The error estimate for the frequenci es (listed in Appendix) is 
obtai ned from the expression given by Lucy & Sweenevl d!971l) 
andlMontgom erv & O'Donoghud d!999l) : 



e(f) 



V6 o- 1 



(3) 



where <x is the standard deviation of the residuals, A expresses 
the semi-amplitude associated to the frequency, AT represents 
the observation span of time and No, the number of points of the 
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binned sample (684). However, these form al errors are often un - 
derestimated compared to the actual ones (Degro ote et al.l2 009). 
The list of frequencies given in the Appendix contains 150 val- 
ues and were stopped at the last significant (<0.01) frequency. It 
is noted that the first 19 frequencies are highly significant against 
both null hypothesis (white and red noise) but also present in 
both partitional datasets. 



Table 1. Table of main frequencies 
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Fig. 6. Evolution of the relative variance of the prewhitened data 
as a function of the number of frequencies already detected by 
the multiperiodic algorithm. 



4. Discussion 

4.1. The inclination of the orbit of Plaskett's Star 

Even though the variations between the maximum and minimum 
observed on the entire run are close to ~ 0.1 mag, the unprece- 
dented quality of the CoRoT light curve gives us enough infor- 
mation to attempt a study of the inclination of Plaskett's Star. 
For this purpose, we use the NIGHTFALL prograrrQ to fit the vari- 
ations due to the orbital motion. This program is based on a gen- 
eralised Wilson-Devinney method assuming a standard Roche 
geometry. 

We have used the information from our Fourier analysis by 
combining the amplitudes and phases of the peaks detected at 
the orbital frequency and its first two harmonics (/ = 0.139 
and / = 0.208 d _1 ). The photometric variations of HD 47129 
tied to the orbital cycle have a peak-to-peak amplitude of about 
19mmag. The orbital light curve displays a broad minimum 
roughly centered on phase 0.0 (the conjunction with the pri- 
mary star being in frontfl This is followed by a broad maximum 



1 For details, see the NIGHTFALL User Manual bv lWichmanH d!998t) 
available at the URL: http : //www.hs .uni-hambourg . de/DE/Ins/ 
Per/Wichmann/Nightfall .html 

2 With respect to the time of phase 0.0 defined in Linder et al. ( 2008), 
the minimum of the light curve is shifted by about 0.03 in phase. This is 
somewhat larger than what can be accounted for by the uncertainty on 
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around phase 0.3, whilst there is no clear secondary minimum, 
but rather a kind of plateau before the brightness decreases to- 
wards phase 0.0. These variations are clearly not due to grazing 
eclipses and cannot be explained by the sole effect of ellipsoidal 
variations either. Indeed, pure ellipsoidal variations in a system 
with a circular orbit, such as Plaskett's Star, would produce a 
light curve with two equally deep minima centered on phases 0.0 
and 0.5. We note that ellipsoidal variations in an eccentric binary 
could potentially acc ount for the obser ved light curve (see e.g. 
the case of Tr 16-1 12, lRauw et al.l2009l). However , the radial ve- 
locity curves of Plaskett's Star (Lin der et al .120081 and references 
therein) provide no evidence whatsoever for a non-zero eccen- 
tricity. An alternative possibility to account for the shape of the 
orbital light curve of Plaskett's Star is to assume a configuration 
where one of the stars has a hot spot on its surface. Indeed, the 
presence of such a hot, and hence bright, region could counter- 



the orbital period quoted by Linder et al. Here we have chosen to take 
phase 0.0 at the time of the minimum of the light curve. 
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balance the small ellipsoidal variations. We fit the light curve 

r~l 1 | I I I | | I I 1 I | I rz 

6.04 h 

6.07 E- 

6.08 ^ 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 

0.5 1 

<p (P - 14.396257 days) 

Fig. 7. Example of a model fit of the orbital part of the light curve 
of Plaskett's Star. The parameters of the model (red thin curve) 
are i = 67°, fill P = 0.54, fills = 0.450, longitude of the spot = 
350°, radius of the spot = 7° and temperature ratio (spot/primary 
star) =1.97. 

of Plaskett's Star with the following assumptions: the mass ra- 
tio as well as t he temperatu r es of the stars are fixed according 
to the results of Linder et al. (2008). A priori, the hot spot could 
be located on the surface of either star. However, to fit the ob- 
served light curve, a spot on the secondary, would have to be 
on the rear side of the star (i.e. turned away from the primary), 
whilst a spot on the primary star would be located on the side 
of the star facing the secondary. The latter configuration seems 
more plausible, especially in com parison with the wind inter- 
action model of Li nder et ID ((2008). We therefore assume that 
the hot spot lies on the surface of the primary. The secondary is 
taken to rotate asynchronously with a rotation period one quarter 
of the orbital period, although we stress that this assumption has 
little impact on our results. At first, the inclination, the Roche 
lobe filling factors of both stars, the longitude of the hot spot, its 
radius and temperature ratio are taken as free parameters. The 
small amplitude of the light curve prevents us, however, from 
deriving strong constraints on the values of all these parameters. 
This dilemma can be illustrated by considering for instance the 
orbital inclination. As a matter of fact, when systematically ex- 
ploring parameter space, we find equally good fits for values of 
i ranging from 30 to 80°. Obviously to reproduce the observed 
light curve with i = 80° requires extremely low filling factors 
of at least one of the stars and this is unlikely since the spec- 
troscopic investigation revealed a luminosity ratio of about 1.9 
for st ars with essentially identical temperatures (see Linde r et"aT] 
2008). Conversely, the lower inclination solutions yield unreal- 
istically large masses. When we constrain the filling factors in 
such a way as to reproduce the spectroscopic luminosity ratio, 
we find the best fit for an inclination of about 67°. Whilst this 
number seems quite 'reasonable', we stress that it must be taken 
with caution. We also attempt to fit the ellipsoidal variations in- 
dividually (i.e., limiting ourselves to the double of the orbital 
frequency and not accounting for a spot) but, once again, no ex- 
act determination of the inclination has been possible. Indeed, 
the fit of this parameter is directly dependent on the filling factor 
of both stars. 

For all the solutions that we find, the longitude of the bright 
spot is found between 345 and 355° (0° corresponding to the di- 
rection of the secondary star and 90° indicating the direction of 
the motion of the primary). Furthermore, this spot has a temper- 
ature ratio of about 1 .9 with respect to the surface temperature 
of the primary. The radius of the spot is found to be about 8 - 
10° for most solutions, except for those with very high orbital 
inclinations and low filling factors. 



A major problem concerning Plaskett's Star remains the dis- 
tance at which the binary system is situated. Indeed, whatever 
the configuration used to fit the light curve, the parameters used 
by NIGHTFALL indicate that the star should be located between 
2.0 and 2.2 kpc. The increase of the distance would imply that 
Plaskett's Star does not b elong to the Mon OB2 association, as 
was already suspected bv lLinder et al.l (120081) . These authors in- 
deed found a discrepancy between the spectral types of both stars 
and their dynamical masses. If we compare the syste mic velocity 
of Pla skett's Star, estimated at 30.6 ±1.8 kms 1 bv lLinder et al.l 
(2008), with the radial v elocities of other O-type stars situated 
in Mon OB2 association (Ma hy et al.l20 09). it appears that these 
values agree between each other and that Plaskett's Star would 
belong to this association. However, this assumption will only 
be checked with measurements from the future Gaia mission. 

In summary, whilst the light curve does not allow us to estab- 
lish the inclination of the orbit, we find that its shape is consistent 
with moderate ellipsoidal variations altered by the presence of a 
rather bright spot on the primary star, facing the secondary. This 
spot is very probably related to the wind interaction between the 
two stars and could actually be due to shock-heated material that 
is cooling whilst it flows away from the stagnation point. 

4.2. Structures present in the light curve 

The multiperiodic analysis of the CoRoT light curve allowed us 
to detect a significant structure composed of a fundamental fre- 
quency (/ = 0.823 d _I ) and its six harmonics but the nature of 
this signal is not yet established. 

In a binary system such as Plaskett's Star, a possible cause 
of spectroscopic or photometric variability could be the asyn- 
chronous rotation. Indeed, a binary system is in synchronous 
rotation when the angular rotation velocity to and the angu- 
lar velocity of orbital motion Q. are equal. In HD 47129, this 
is clearly not the case: adopting the projec t ed rot ational ve- 
locities and stellar radii from Linder et al ] (120081) . we esti- 
mate rotational periods of about 10.3 (0.1 d _I ) and 1.7 days 
(0.6 d -1 ) for the primary and secondary, respectively. As a re- 
sult of asynchrono us rotation, tidal interactions may create non- 
radial oscillations jj Willems & Aertsll200l iMoreno etail 120051: 
Koenigsberger et al.ll2010l) . These oscillations produce a pattern 
of azimuthal velocity perturbations superposed on the unper- 
turbed stellar rotation field. These tidal flows lead to the dissi- 
pation of energy due to the viscous shear, thereby impacting on 
the surface temperature (and hence brightness) distribution. The 
azimuthal velocity components of tidal interactions are expected 
to produce strong line profile variability that resembles the typ- 
ical signature of non-radial pulsations (bumps in the line profile 
that migrate from the blue wing of the spectral lines to the red 
wing). The most spectacular effects are expected in eccentric bi- 
nary systems and this model was successfully used to explain the 
line profile variabili t y of the eccentric B-type b inary a Virginis 
(IMoreno et alj|2005l: [koenigsberger et al.ll2010l) . Simulations of 
the influe nce of the tidal effects on the radial velocities were 
shown bv I Willems & Aertsl (l2002l) . 

In a binary system with a circular orbit, the non-synchronous 
rotation induces variations of the radius of the non-synchronous 
rotating component with a super-orbital period that is longer 
than both the orbital and the rotational period (IMoreno et al.l 
2005). Whilst these super-periods could leave a photometric sig- 
nature, they should not be visible in the line profile variabil- 
ity whi ch is mostly dominate d by the azimuthal velocity com- 
ponent. IMoreno eta D d2005l) found that, in the circular non- 
synchronous case, the line profiles should only display a few, 
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rather broad bumps. The skewness of the line profile and its ra- 
dial velocity (compared to the unperturbed profile) should dis- 
play two maxima and minima per orbital cycle. The correspond- 
ing frequency should thus be twice the orbital frequency which 
is different from / = 0.823 d" 1 . 

We then consider that the structure composed of 7 frequen- 
cies could be due to the rotation of one component of the binary 
system. In thi s context, we clearly see by adopting the rotational 
periods from iLinder et al.l (|2008) that this set of frequencies is 
unlikely to correspond to the rotation of one star. A similar con- 
clusion is expected if we t ake into account th e projected rota- 
tional velocity evaluated bv lLinder et al.l d20 08) for both stars at 
about 75 and 310 kms" 1 fo r the primary and seco ndary com- 
ponent, respectively. Indeed, Ba gnuolo et al.l d 19921) derived the 
radii from the luminosities and the effective temperatures of both 
stars. They found values close to 14.7 R G and 9.7 R e for the 
primary and secondary stars, respectively. We estimate a rota- 
tional period of about 9.9 and 1.6 days for the primary and sec- 
ondary stars, respectively, corresponding to rotational frequen- 
cies of about 0.1 and 0.6 d . Alternati vely, if we adopt the t he- 
oretical values of the radii quoted by iMartins et alJ (120051) as 
a function of the spectral type of each component, we obtain 

21.1 R Q and 14.2 R Q for the primary and secondary stars, re- 
spectively, by assuming that the system is composed of an 08 1 
and an 07.5 III star. These values yield a rotational period of 

14.2 days for the primary and 2.4 days for the secondary. These 
results suggest that the primary could be in synchronous rotation 
while this is clearly not the case for the secondary, which would 
have a rotational frequency of about 0.42 d~'. In both cases, it 
thus appears that the fundamental frequency found in the Fourier 
analysis is not representative of a rotational motion, except if we 
could invoke a spotted surface for the secondary. But, in this 
case, the lack of odd multiple of the rotation frequency is sur- 
prising. By assuming a rotational frequency of / = 0.823 d _1 , 
we would obtain a radius of about 7.2 R for the secondary star, 
which seems unrealistic for such stars. However, it is possible 
that the frequency at / = 0.368 d" 1 or / = 0.399 d" 1 is repre- 
sentative of the rotational period of the secondary component. In 
the latter case, the equatorial radius of the secondary star could 
be estimated close to 16.7 or 15.4 R & , respectively. 

Next, we suppose tha t this group of freq uencies could arise 
from wind interactions. Wig gs & Giesl (1 19921) reported on vari- 
ations of the strength of the high-velocity wings of the Ha and 
He i A 6678 emission line profiles. The total equivalent width of 
the blue and red wings of Ha was found to display a roughly 
recurrent modulation with a time scale o f 2.82 days (0.35 d" 1 ) , 
although the sampling of the spectra of IWiggs & Giesl {l992) 
is admittedly not sufficient to characterize such rapid varia- 
tions. These authors attributed the variations to instabilities of 
a radiatively cooling wind-wind interaction. It is interesting to 
note that this frequency is very close to the "broad peak" de- 
tected in the analysis of our photometric data (Sect. 13. U . It 
has to be stressed th ough that the XMM-Newton observations of 
ILinder et alJ {2006) did not reveal any indication of X-ray vari- 
ability that could be related to the existence of such instabilities 
in the wind interaction zone. 

Finally, we consider that these frequencies could be gener- 
ated by non-radial pulsations. Indeed, massive stars are com- 
posed of a convective core and a radiative e nvelope, and grav - 
ity as well as acoustic modes can be excited dAerts et alJlioToh . 
Generally, for massive stars situated on the main-sequence band, 
the p-modes propagate with periods of a few hours while the 
g-modes have longer peri ods, of orde r of da ys. However, in 
the case of Plaskett's Star. ILinder et al.l d2008l) concluded at an 



evolved system in a post Roche lobe overflow evolutionary stage. 
The inner structure of such stars is different from main-sequence 
stars, which decreases the g-modes periods and, at the same 
time, increases the p-modes perio ds (|Handlerll2005l) . Another 
family of modes, the strange modes dSaio et al.l 1998b . could also 
play a role in the interpretation of this structure of frequencies. 
These modes are known to have a propagation zone near the sur- 
face of massive stars. The exact determination of these modes 
requires theoretical models to predict variations produced by 
non-radial pulsations in a rapidly rotating massive star as well 
as an intense spectroscopic monitoring to characterize the vari- 
ability of the line profiles. However, for Plaskett's Star, massive 
main- sequence star models computed with ATON dVentura et al.l 
2008) and an initial metallicity of 0.015, a range of mass be- 
tween 30 and 70 M Q and different mass-loss rates show that 
non-radial pulsations with frequency of the order of 0.8 d _1 can 
indeed be found. This frequency could be generated by modes 
with I = 2, 3 or 4. 

In summary, we conclude that the alternative scenarios 
linked to rotation or to winds, in both stars, fail to explain the 
structure of the Fourier spectrum, especially the dominant fre- 
quency and its set of harmonics. The assumption of multiperi- 
odic non-radial pulsations remains the most plausible one to un- 
derstand the low-frequency variability detected in the CoRoT 
light curve of Plaskett's Star. 

5. Conclusions 

We presented the photometric analysis of the very massive bi- 
nary system, HD47129, observed by CoRoT during the second 
run SRa02 (~ 34.33 days). The light curve shows indications 
of intrinsic variations of the stars, providing evidence of a large 
number of frequencies. We extracted a total of about 43 frequen- 
cies (TableQ}, significant under a null hypothesis of red noise, by 
using two different techniques based on a standard prewhitening 
and a multiperiodic algorithm. We emphasize that all the fre- 
quencies, reported in the present paper, have not necessarily a 
physical meaning. 

This analysis highlighted a group consisting of one funda- 
mental frequency (0.823 d _1 ) and its six harmonics. In addition, 
a second structure formed by the orbital frequency (0.069 d -1 ) 
and its two harmonics allowed us to investigate the light varia- 
tion due to orbital effects. This analysis revealed the presence of 
a hot spot, probably located near the primary star, facing the sec- 
ondary component and having an origin in the colliding wind in- 
teraction zone. However, the determination of the exact value of 
the system inclination turned out to be ambiguous. A third struc- 
ture, formed by a much wider peak and thus certainly composed 
of two frequencies at 0.368 d _I and 0.399 d _1 , has also been 
detected. One of these frequencies could be related to th e rota- 
tion of the secondary even though Wigg s & Giesl d!992h rather 
suggested an origin in the winds of the stars. 

The future work will be devoted to the interpretation of these 
frequencies in terms of asteroseismology in order to remove an 
additional part of mystery concerning Plaskett's Star. Moreover, 
a study of the line profiles either of the secondary or of the 
primary component, obtained from an intense high-resolution 
spectroscopic campaign, could provide further constraints on the 
properties of the pulsation modes. Indeed, the broad and shallow 
line profiles of the secondary star favor this detection although 
the orbital motion and the presence of the primary star will cer- 
tainly render this task more difficult. 
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Table .2. Total list of frequencies. The errors given on the frequencies, semi-amplitudes and phase correspond to the l-<x error. The To of the phase 
is equal to 2,450,000. The last column marks the frequencies present both in the first part of the light curve and in the second one. 
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